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Abstract We investigate the influence of large-scale meridional circulation on so- 
lar p-modes by quasi-degenerate perturbation theory, as proposed by Lavely and 
Ritzwoller (1992). As an input flow we use various models of stationary meridional 
circulation obeying the continuity equation. This flow perturbs the eigenmodes of an 
equilibrium model of the Sun. We derive the signatures of the meridional circulation 
in the frequency multiplets of solar p-modes. In most cases the meridional circula- 
tion leads to negative average frequency shifts of the multiplets. Further possible 
observable effects are briefly discussed. 
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1. Introduction 

Meridional circulation is a large-scale flow observed on both hemispheres of the 
solar surface (Duvall, 1979; Hathaway, 1996; Komm, Howard, and Harvey, 1993). 
Its predominant direction is from the equator to the poles, and its amplitude is of 
the order of 15m/s. As mass does not accumulate in the polar regions a return 
flow from the poles to the equator is suspected deeper within the solar interior. The 
top half of the convection zone contains approximately 0.25% of the solar mass, the 
mass of the bottom half is approximately five times larger. Consequently, a poleward 
flow of lOm/s in the top half of the convection zone could be compensated by an 
equatorward flow of 2m/s in the lower half. The transport of magnetic flux from 
mid to low latitudes by such a flow at the bottom of the convection zone would last 
approximately 10 years, which is close to the period of the solar magnetic cycle. 

In addition to magnetic flux, the meridional flow also transports angular mo- 
mentum. Indeed, the circulation played a key role in early theories of the solar 
non-uniform rotation as well as of the magnetic cycle (Bjerknes, 1926; Kippenhahn, 
1963). More recently, differential rotation has been explained in mean- field models as 
a consequence of the Reynolds stresses (Riidiger, 1980; Kiiker and Stix, 2001; Kiiker 
and Riidiger, 2005), and in three-dimensional numerical models by the influence 
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of the Coriohs force on global convection (Miesch et a/., 2000; Miesch, Brun, and 
Toomre, 2006). Nevertheless, meridional circulation occurs as well in these models, 
and in solar-cycle models it has regained popularity, since the mean-field latitude 
migration along the surfaces of isorotation that occurs in the traditional aQ dy- 
namo (Parker, 1955) does not seem to suffice. The effect of the circulation on the 
butterfly diagram had been demonstrated by Roberts and Stix (1972); it is consid- 
ered to be essential in more recent versions of the aft dynamo, which therefore have 
been termed 'flux- transport dynamos' (Choudhuri, Schiissler, and Dikpati, 1995; Dik- 
pati and Charbonneau, 1999; Nandy and Choudhuri, 2002; Rempel, 2006; Rempel, 
2006). 

Local helioseismology has investigated the strength of the meridional flow in the 
solar interior. By means of ring diagram analysis (Hill, 1988) Haber et a/., (2002) 
inverted data for the circulation in a 15 Mm deep region below the solar surface. In 
data from 1998 they found a flow emerging at high northern latitudes with equa- 
torward orientation. This was interpreted as an evolving second cell of circulation. 
Further studies on the evolution of the flow either by ring-diagrams and by time- 
distance helioseismology (e.g. Zhao and Kosovichev, 2004, Zaatri et a/., , 2006) show 
predominantly a polward flow with a strong variability in the outer 15 Mm of the 
Sun. The velocity reaches 40m/s. These flndings were interpreted as the upper parts 
of meridional circulation cells in the two hemispheres. 

Theoretically, the influence of global-scale stationary flows on solar p-modes was 
studied in detail by quasi-degenerate perturbation theory (Lavely and Ritzwoller, 
1992). These studies were successfully used to solve the forward problem for the 
influence of differential rotation on the p-modes. The results lead to an improved 
inversion method for determining the radial dependence of the differential rota- 
tion (Ritzwoller and Lavely, 1991). Following Lavely and Ritzwoller (1992), Roth 
and Stix (1999) studied the influence of large-scale sectoral poloidal flow compo- 
nents that could be related with giant convection cells. They found that such flows 
yield additional frequency shifts that can only be described with quasi-degenerate 
perturbation theory, as these frequency shifts are effects of higher order. In a sequel 
study Roth, Howe, and Komm (2002) were able to show that sectoral poloidal flows 
could only be found with the current inversion methods of global helioseismology 
as long as they exceed an amplitude of lOm/s. The meridional flow was found to 
be not detectable by the current inversion methods as the frequency splittings are 
fltted by an incomplete set of basis functions, which are tailored to measure only zonal 
toroidal flows. However, no detailed study on the effect of the meridional circulation 
on the oscillation frequencies was given, and few other attempts to derive observable 
signatures of the meridional flow in global helioseismology data exist (Woodard, 
2000). 

In this contribution we concentrate on a theoretical study of the influence of 
the meridional circulation on the solar p-mode frequencies and describe a possible 
observable effect. This effect is signiflcantly smaller than the frequency splitting 
caused by solar differential rotation. But as time series of global oscillation data 
exist that cover more than 10 years, the necessary frequency resolution might be 
available to detect it. The advantage of studying the meridional circulation by global 
helioseismic techniques is a possible inference of information from greater depths. 
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2. Frequency shifts caused by meridional circulation 

The effect of the meridional circulation on solar oscillations shall be investigated 
by solving the forward problem. We use quasi-degenerate perturbation theory as 
proposed by Lavely and Ritzwoller (1992) to calculate shifts of the oscillation fre- 
quencies. A short outline of the mathematics is given in the follwing. 

As described by Lavely and Ritzwoller (1992) and Roth and Stix (1999) we con- 
sider a spherically symmetric equilibrium model of the Sun. For this purpose "Model 
S" from Christensen-Dalsgaard et a/., (1996) is used. The solar p-modes are adiabatic 
eigen-oscillations ^/^ of small amplitude, where k stands for the three indices har- 
monic degree /, azimuthal order m and radial order n. We calculate the eigenfrequen- 
cies and eigenmodes numerically with the ADIPACK code (Christensen-Dalsgaard 
and Berthomieu, 1991). 

The equation governing the eigenmodes is 

^oik = -Po^fc^/c ^ (1) 

where po is the density and uo^ the oscillation frequency The operator Cq acts on 
the eigenoscillation by 

£0^ = - VP' + pQg' + pgo , (2) 

where the primed quantities are the Eulerian variations of pressure, gravitational 
acceleration, and density caused by an oscillation mode (e.g., Stix, 2004, Unno et al., 
, 1989). 

Describing the disturbing effect of the meridional circulation we replace in Eq. (1) 
the operator Cq^ the squared eigenfrequency cj^, and the eigenfunction with 

Co Co -\- Ci , 

^ ojj , (3) 

~^ ' 

where the perturbation operator Ci is defined in terms of the meridional circulation 
uo and acts on the eigenmodes as 

Cii^k) = -^ioUrefPoiuo • V)^fc . (4) 

The perturbed eigenfunctions are expressed in terms of the unperturbed normal 
modes of the equilibrium model 

keK, 

where the index k and the subscript j are elements of the subspace JC that is spanned 
by the eigenfunctions which are quasi- degenerate. We identify j with that k for 
which the expansion coefficients are maximal in magnitude. Following Lavely and 
Ritzwoller (1992) we use 

= ^vei + A (6) 

for the perturbed frequency The reference frequency cjref should be chosen in the 
vicinity of the eigenfrequencies of the modes in the subspace JC. 
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According to Lavely and Ritzwoller (1992) the problem can be turned into an al- 
gebraic eigenvalue problem by substituting (5) and (6) into the perturbed eigenvalue 
problem, and making use of the orthogonality of the eigenfunctions 

aiZk'k = ^ a^k^jh'k for k' e IC . (7) 

keic keJC 

The elements of the matrix Z are given by 

Zk'k l^^n^'/ i^rei ^k)^k^k'k^ , (8) 

with 

^;?l> = -y'^fc'-A4^'r. (9) 
The normalization of the eigenfunctions is given by 

Nk = I ^ PoKr + /(/ + dr . (10) 

2.1. The model of the meridional circulation 

We use a linear superposition of Legendre polynomials Ps (cos 6) with degree s for 
our representation for the meridional velocity field 

ito = ^^^Us{r)csPs{cos9)er -\- Vs(r)\/ hCsPs (cos 9) , (11) 

s 

where Vh is a horizontal gradient and is the colatitude. The factor Cs is used for 
normalization of the Legendre polynomials and is determind by = (2s+l)/47r. The 
expansion coefficients Us{r) and Vs{r) are functions of the radial coordinate r, and 
represent the radial and horizontal strength of the fiow component with degree s. In 
this model the meridional velocity field is a zonal poloidal fiow which is independent 
of longitude. 

The meridional circulation shall be free of divergence, i.e. V • (pouo) = 0, which 
enables to express Vs in terms of Us, 

pors{s + l)vs = dr{r^pQUs) . (12) 

We adopt a simple model for the depth dependence of Us{r)^ which is given by 

Us{r) = Asin(nc7r— ) for < r < Rq , 

Us(r) = otherwise . (13) 

The return fiow of the meridional fiow is closed at the radius given by the parameter 
r5, which we set equal to the bottom of the convection zone = 0.7137?© (Basu and 
Antia, 1997). The number of circulation cells in depth is defined by the parameter 
Tie. We select the amplitude A such that the horizontal fiow component Vs at the 
solar surface has a maximum amplitude of 15 m/s. 
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Below we present results of the effects on the solar p-mode frequencies of 15 various 
models of the meridional circulation. We vary s between 2 and 6 and ric between 1 and 
3. The higher degrees s represent more than one cell in each hemisphere; such multi- 
cell circulation occurs in numerical simulations (Brun and Toomre, 2002; Miesch, 
Brun, and Toomre, 2006), and has been observed by means of sunspot tracing on 
the solar surface (Tuominen, Tuominen, and Kyrolainen, 1983; Wohl and Brajsa, 
2001). The models with odd s have an equator-crossing flow; we include these as 
such asymmetric flows have been found on the Sun as well (Haber et al., 2000; Zhao 
and Kosovichev, 2004). Figure 1 illustrates the 15 models. 

Due to Eq. (12) models with higher degree s show stronger radial flows in the 
down- and upfiow channels whereas the horizontal flow strength is the same in all 
models. 



2.2. Theoretical frequency shifts 



We evaluate the integration of (9) by using the model of the meridional circulation 
(11) and (12). The integral over longitude does not vanish only if = m. The result 
is 



l{m) = 8icJref7r(- 1)^^7/7/' 

Rq 

Us Rs — dr 





Hs 



s(s-hl) 



por dr 



s I I' 
m —m 



1(14) 



The integral kernels Rs (r) and Hs (r) are the equivalents of the more general poloidal 
flow kernels given in Lavely and Ritzwoller (1992) 



Rs{r) = 



2 - ^^'^ ) + (,0 1-1 

[Z(Z + l)-Z'(Z' + l)](l + (-l)^+'+'') 



-5;.^.(i + (-i)^+'+'')n;.n. (^^ [ 



S I I 

1 -1 



(15) 



where = \/ (2x -|- l)/47r and Vix = a/x(x -|- 1)/2. The 3x2 arrays are Wigner-3j 
symbols that result from the coupling of the angular momenta of the oscillations 
and the flow involved (Edmonds, 1974). Due to this coupling of angular momenta 
the matrix elment H^'n^l'lij^) is non- vanishing if certain selection rules are fulfilled. 
The first rule arises from properties of the Wigner-3j symbols which vanish exept 
when the harmonic degrees /, and s satisfy a triangular condition. The second 
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Figure 1. Streamlines for the flows puo of the 15 meridional circulation models as functions of 
solar colatitude and radius. Black indicates counterclockwise flow, grey clockwise. The arrows 
indicate the flow velocity at some positions. The number of cells in horizontal direction changes 
from top to bottom between 2 and 6. The number of cells in depth changes from left to right 
between 1 and 3. 



rule follows from (15); the sum of the degrees must be even otherwise the factor 
(1 + (-1)^+^+^') in Rs and Hs vanishes, 



s + / + = mod 2 . 



(16) 
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We can simplify Eq. (14) by defining a poloidal flow kernel Ks{r) 

Rq 

^n'n,Z'zM =: 8icJref7r(-l)'^7z7z' X^75 J UsKs{r)por'^ dr (^^^ ^ ^ • (^^) 

^ 

According to this result the meridional circulation leads to an effective shift of 
the eigenvalues, i.e. the mode frequencies, if the selection rules are fulfilled. These 
shifts are given by the eigenvalues A of the matrix Z, or in more detail, the corrected 
mode frequency Cunlm is given by 

(18) 

Following the concept of quasi-degenerate perturbation theory, the coupling of only 
two modes in the presence of only one flow component was investigated by Roth and 
Stix (1999). They found that due to couphng of two modes by any poloidal velocity 
field the absolute magnitudes of the frequency shifts are equal for both couplers, 
but the sign is negative for the mode with the lower frequency. In the case of two 
multiplets where pairs of modes are coupling the frequency shifts are a function of 
the azimuthal order m. The functional dependence on m is given by the absolute 
magnitude of the Wigner-3j symbol as a function of m. It then follows that the 
shifts within one multiplet are either always negative or always positive. From this 
it follows for the case of the meridional circulation that due to the dependence of 
the Wigner-3j symbol on m the frequency shifts are always symmetric to the mode 
with m = which in general is also shifted. 

We would also like to point out that due to the second selection rule (16) the flow 
kernel Ks vanishes if l' = /, i.e. there are no contributions to the diagonal entries 
on the matrix Z from the matrix elements Hn'n,in{'^)- As the matrix elements 
Hn'n,l'l{'^) are complex and Z is Hermitian, a change of the orientation of uq does 
not change the frequency shifts. 



3. Results 

The result of our investigation are frequency shifts for the mode multiplets with 
< / < 300 and 1 < n < 30 calculated numerically by setting up the full coupling 
matrix (8), i.e. taking all possible couplings of modes into account, and evaluating 
the eigenvalues. 

We first focus on the frequency shifts in the multiplets. According to the properties 
of the Wigner-3j symbols the shifts within a mulitplet are symmetric to the mode 
with m = 0. Figure 2 displays an example. The frequency splitting is stronger for 
meridional circulation components with higher degree s. In the example of Fig. 2 
the shifts are in the order of 1 /xHz for s = 6. The overall effect of the superposition 
of the meridional circulation components results in shifts in the order of 3 /xHz away 
from the unperturbed frequency. The difference between the summed-up effect of 
models with nc = 1 cell in depth and models with nc = 3 cells in depth is in the 
order of 0.1 /xHz. 

To give a better overview of the effect on all modes, Fig. 3 displays the mean 
frequency shifts of the multiplets as a function of the unperturbed mode frequencies. 
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Figure 2. Splittings of the p-mode multiplet n = 1, I = 292, u = 2263.08 /xHz as a function of 
azimuthal order m. The splittings are caused by different meridional flow models: s = 2, Uc = 1 
(purple), s = 3, Uc = 1 (blue), s = 4, Uc = 1 (green), s = 4, ric = 1 (yellow), s = 4, ric = 1 
(red). The total effect of a sum of these flve flow models is displayed by the thick red line, the 
black line on top of this gives the result for a sum of the flve flow models but with Uc = S. 
The unpertubed mode frequency is given by the dashed line. 

The frequency shifts were obtained for the single meridional flow components. One 
common result is a negative frequency shift of most modes. We find positive frequency 
shifts only in a few multiplets. The origin of this effect is given by the shape of the l-v 
diagram. According to quasi-degenerate perturbation theory the multiplet nearest 
in frequency fulfilling the selection rules causes the strongest shift (Roth and Stix, 
1999). As shown for the case s = 2 in Fig. 4 in most cases these nearest neighboring 
modes have a higher frequency. This is explained by the first selection rule and the 
curvature of the ridges in the l-iy diagram. Because of the low degrees s of the fiow 
components, the difference of the harmonic degrees of two couphng modes |/ — l^\ 
can not be larger than s. Possible coupling partners for a particular mode can then 
only come from a very narrow region in /. As the slopes of the ridges are steeper 
towards lower harmonic degrees, the frequency spacing between the modes decreases 
along a ridge. In addition, the frequency difference between two ridges increases 
with /. Therefore the nearest neigboring mode fulfilling the selection rules usually 
lies on the same ridge and has a higher frequency and a higher harmonic degree. 
Figure 5 shows that in most of the investigated cases the coupling partners lie on 
the same ridge. Interestingly there is some structure visible in these plots. There are 
areas where all coupling partners are coming from the same ridge or from one ridge 
above or below. These areas are defined by the curvature of the ridges. With higher 
degree s of the meridional fiow the coupling of modes from different ridges becomes 
more frequent. These plots indicate only the position of the nearest coupling partner, 
because this partner causes the strongest effect. The effects of coupling partners with 
greater distances in frequency are weaker. 

Another common result is that the frequency shifts are strongest for the modes 
with low radial order n and with high harmonic degree /. This is again explained 
by the location of the nearest neighbor in frequncy (cf. Fig. 4). The lower the radial 
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Figure 3. The mean frequency shifts of the multiplets caused by the 15 different configurations 
of the meridional circulation (cf. Fig. 1) as functions of the unpertubed mode frequencies. The 
color code gives the radial order n of the modes. The degree of the meridional circulation 
affecting the modes changes from s = 2 (top) to s = 6 (bottom), the order Uc of the meridional 
circulation changes from 1 (left) to 3 (right). 



order n the smaller the frequency difference to the next mode. Along a ridge the 
frequency difference to the next mode decreases causing a stronger splitting of the 
multiplets. On the ridges with high radial order, the frequency spacing between two 
neigbouring modes is larger than on the ridges with low radial order. 

These two common results apply in principle to all poloidal flows, not only to 
the meridional circulation. Therefore we can conclude that the overall effect of large- 
scale poloidal flows in the convection zone is a lowering of the mode frequencies. This 
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Figure 4. The frequency difference to the next nearest mode fulfiUing the selection rule for 
the meridional circulation component with s = 2. The frequency difference is plotted versus 
mode frequency. The color code gives the radial order of the modes. 




Figure 5. Theoretical l-u diagrams with colors indicating the origin of the coupling partner 
with nearest frequency. Green: coupling partner from the same ridge; blue: coupling partner 
from one ridge above; red: coupling partner from one ridge below, black: coupling partner is 
more than one ridge apart. The degrees s of the meridional circulation components change 
from 2 to 6. 
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effect could contribute to reduce the discrepancy between the theoretically calculated 
mode frequencies and the observed frequencies. A similar result was already found 
by (Zhugzhda and Stix, 1994; Stix and Zhugzhda, 1998). 

Investigating the results in detail for the used meridional circulation models we 
find that the frequency shifts are increasing from s = 2 to s = 6. In the case 
of s = 2, nc = 1 the total average of theses shifts for all modes affected by the 
meridional circulation is 0.01 /xHz, whereas in the case s = 6, ric = 1 the magnitude 
of the average shift is 0.1 /jHz. 

There is also a small difference noticeable in the frequency shift if the number 
Tie of cells changes. The difference of the average shifts calculated for the models 
s = 2, nc = 1 and nc = 3 is -0.002 /xHz. Whereas the difference of the average shifts 
determined for the models s = 6, nc = 1 and nc = 3 is -0.03 /iHz. However, the 
scatter in the frequency shifts is high. In the case of s = 2, nc = 1 the largest shift 
observed is 22.5 /xHz, and in the case of n the case of s = 6, nc = 1 the largest 
shift observed is 79.4 /xHz. Strong shifts above 10 /xHz occur for the modes with high 
harmonic degree /. 

4. Discussion 

In this paper we used simple models of the meridional circulation to investigate their 
influences on the solar p-mode frequencies. The simplest model consisted of one cell 
per hemisphere and depth with a maximum horizontal flow velocity of 15m/s on the 
solar surface. The most complicated model we used had three cells per hemisphere 
and three revolution cells in depth with a horizontal flow velocity of 15m/s at 
the surface, too. We performed numerical calculations based on quasi-degenerate 
perturbation theory to obtain the frequency splittings of the solar p-modes due to 
this meridional circulation models. We find that the meridional circulation lifts the 
degeneracies of the multiplets. For the simplest model the shifts are on average only 
0.01 /jHz with a few shifts up to several /iHz. Models with more cells per hemisphere 
affect the p-modes stronger. For the model with three cells per hemisphere we find 
an average shift of 0.1 /iHz, with many shifts in the order of 1 /xHz. In most cases 
the shifts are negative due to the fact that the next neigboring mode in frequency 
dominates the shifting; and due to structure of the l-v diagram this nearest neighbor 
has usually a higher frequency causing therefore negative shifts. 

Comparing models with the same number of cells in latitude but different number 
of cells in depth we find only tiny differences in the resulting frequency splittings. 
In the case of s = 6 the difference in the frequency splitting is on average 0.03 /xHz. 
However a number of differences in the order of 1 /iHz can occur. 

On the Sun the meridional circulation consists probably of a superposition of flow 
components with different numbers of cells in depth and latitude. In contrast to the 
models used in our work, the amplitudes of these flow are in reality likely to be very 
different and highly variable in time. Nevertheless, based on our results we are opti- 
mistic that the meridional circulation on the Sun could leave an observable signature 
in the p-mode frequencies. In order to detect the effect a frequency resolution of at 
least 0.1 fiHz must be achieved. This could be done by a several month long time 
series. In order to be able to distinguish not only the various components of the 
meridional circulation in the orders s but also in the radial orders nc, several years 
of data need to be averaged to obtain the required precision in the frequencies. As 
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such long data sets exist it will be worth while to search for this effect in the p-mode 
frequencies. 

Compared to the splitting caused by the differential rotation the effect of merid- 
ional circulation is small. But as the effect of the rotational splitting is odd it can 
be separated from the even meridional splitting. However, e.g. asphericities and the 
magnetic field might cause symmetrical frequency shifts, too. Therefore, a lot of 
forward modelling will be necessary before the effect of the meridional circulation can 
be disentangled from these other effects. In this sense, one possible future extension 
of our work is the determination of the frequency shifts due to more sophisticated 
models of the meridional circulation, e.g. from three-dimensional numerical models. 
This might allow tailoring inversion routines for estimating the meridional circulation 
in the Sun from global solar oscillation frequencies. 
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